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In models with large extra dimensions, black holes may be produced in high-energy particle 
collisions. We revisit the physics of black hole formation in extensive airshowers from ultrahigh- 
energy cosmic rays, focusing on collisional QCD and black hole emissivity effects. New results for 
rotating black holes are presented. Monte Carlo simulations show that QCD effects and black hole 
spin produce no observable signatures in airshowers. These results further confirm that the main 
characteristics of black hole-induced airshowers do not depend on the fine details of micro black hole 
models. 
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I. INTRODUCTION 



In an effort to unify gravity with the other fundamental forces, we are faced with the hierarchy problem. The 
fundamental scale of gravity is about 17 orders of magnitude higher than the TeV scale, where electromagnetic and 
weak forces unify. The hierarchy problem may be solved by the introduction of large extra dimensions [l| . In this 
model, the Planck scale Mp\ is related to the fundamental scale of gravity ~ 1 TeV by the relation Mpj ~ V n M™ +2 , 
where V n is the volume of the extra n-dimensional space. Gravity is a strong force in the higher-dimensional spacetime 
but appears weak to a four-dimensional observer due to its "leakage" in the extra dimensions. Gravitons may propagate 
in all dimensions (bulk). Compatibility with known sub- TeV physics restricts the propagation of all Standard Model 
(SM) fields to three spatial dimensions (brane). 

■ One of the effects of the increased strength of gravity would be the production of TeV-scale Black Holes (BHs) in 
high-energy particle collisions Q. Micro BHs could be produced in man-made particle colliders, e.g. the Large Hadron 
Collider 0, HJaL or naturally in Earth's atmosphere by Ultra High Energy Cosmic Rays (UHECRs) interacting with 
air nucleons [6j, M, 10. (For reviews, see Refs. (lQ|-) Once formed, these BHs would immediately decay through loss 
of excess multipole moments (balding phase), Hawking emission [llj (evaporation phase) and final n-body decay or 
remnant production (Planck phase). SM fields may be originated in each of these stages, providing a means to detect 
the BHs. For atmospheric events, the visible imprint would be an extensive airshower initiated by these SM quanta. 

The characteristics of BH-induced airshowers can be investi gate d with Monte Carlo techniques. A Monte Carlo 
code for BH formation and airshower generation is GROKE |12j. GROKE simulation of BH events proceeds in 
three stages. First, the BH is formed by the collision of an UHE neutrino and an air nucleon parton. Some of the 
Center-of-Mass (CM) energy is lost in the process (about 40% for head on collisions, monotonically increasing with 
the impact parameter pH, Il4j). SM unstable particles from the BH decay and the nucleon remnant are hadronized 
using a high-energy physics program for event generation (PYTHIA) [15]. PYTHIA's output is then injected into 
a simulator of extensive airshowers (AIRES) [16]. Simulations show that BH airshowers generally rise faster, have 
broader peak and higher variation in the total energy than SM airshowers because of the "democratic" nature of BH 
decay. BH events are also characterized by a larger muonic content at the ground compared to SM events due to the 
dominant hadronic channel in the BH evaporation phase. A complete discussion of BH signatures can be found in 
Refs. H,[l2[. 

BH searches require the identification of observational signatures that do not depend on the fine details of the 
model. This can be achieved by improving the theoretical description of the event and then tes ting the stability 
of the airshower characteristics against these theoretical refinements. Previous investigations [1, Il2j ] neglected or 
approximated various aspects of the physics of BH formation and decay such as QCD effects, BH spin and particle 
emissivities. QCD effects in the fragmentation process may lead to changes in the amount of visible energy deposited 
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in the airshower by the nucleoli remnant. Changes in particle emissivities due to spacetime dimensionality or BH 
rotation may affect rapidity, peak variation and muon content of airshowers. The aim of this paper is to check the 
stability of BH airshowers characteristics when these effects are included. The result of our investigation is that 
inclusion of collisional QCD effects, changes in particle emissivity and BH rotation do not significantly affect the BH 
airshower development: Observational signatures of BH events are robust. Natural units are used throughout the 
paper with Af* = 1. 



II. QCD AND EMISSIVITY EFFECTS 



QCD effects in airshower generation and development include initial- and final-state radiation, fragmentation, and 
the hadronization process of nucleon remnant and unstable quanta. Since PYTHIA is designed to handle initial- 
and final-state radiation, multiple scattering, beam remnant and hadronization, these effects can be investigated by 
modifying the GROKE code |12j to include the BH event as a PYTHIA external process. 

The BH airshower is initiated by the decay products of the BH, the nucleon remnant and jets from initial- and 
final-state radiation. The colliding parton is taken from the Parton Distribution Functions at a very high energy 
scale (> TeV). This implies that it typically emits quite hard initial-state radiation before it collides to form a black 
hole, resulting in additional jets. In addition, the nucleon remnant undergoes soft and semi-hard multiple scatterings 
which will contribute to the airshower. Previous investigations of BH airshowers did also not take into account color 
conservation. However, the nucleon remnant is color connected to the BH decay products; the color flow will also 
hadronize into more jets. In our investigation, color flow is implemented in the Nq — > oo limit of QCD [l5j . 

Total multiplicity (N) and multiplicity per species (iVj) of the BH decay are essential to determine the airshower 
characteristics. Earlier studies [1, [l2j used approximated thermally-averaged emissivities (graybody factors) for the 
evaporation phase. Recently, exact graybody factors for higher-dimensional nonrotating BHs were calculated in Ref. 
[l7| . Changes in greybody factors are specially relevant for higher-dimensional spacetimes, where graviton emission 
is highly enhanced. These results are implemented in GROKE following Ref. The total multiplicity in the 
evaporation phase is 

(n + l)S Ejfig^ 

4n V,r/P ; rr • [! 

where S is the initial entropy of the BH, Cj are the degrees of freedom of the z-th species, and T-p i and Tfi i arc 
the fraction of radiated power and the emission rate per degree of freedom, respectively. The decay multiplicity per 
species is 

Ni = N *^ . (2) 

The use of exact greybody factors (nonrotating case) leads to a slight reduction in the output of visible energy and 
an enhancement of graviton multiplicity in the evaporation phase compared to previous studies. These effects are 
generally of order ~ 1 or less. 

QCD and emissivity effects on the BH airshower development can be determined by looking at the longitudinal 
development of the e + e~ component of the airshower and the muonic content at ground. The left panel of Fig. [1] 
compares 50 BH airshowers with and without QCD and emissivity effects (primary neutrino energy E u = 10 19 eV, ten 
spacetime dimensions). The average depth of the airshower maxima X m is not significantly affected by the inclusion 
of initial- and final-state radiation, color conservation and exact emissivities. This can be qualitatively explained by 
looking at the energy distribution of the BH airshower initiatiors after the hadronization and fragmentation process. 
PYTHIA's output shows that the additional jets from initial- and final state radiation are generally too soft to affect 
the airshower development, which is mainly determined from the evolution of the nucleon remnant and the hard 
hadronic jets from the BH evaporation. The implementation of color conservation slightly changes the details of the 
hadronization process of previous studies. However, the main characteristics of the airshower depend on the hadronic 
nature of the event rather than the details of the fragmentation. Differences in the fragmentation model are washed 
out by uncertainties in the airshower development. A similar qualitative explanation applies to the emissivity effects. 
Even for massive BH events, when quanta from the evaporation phase dominate over quanta from the Planck phase 
and the nucleon remnant, changes due to the use of exact greybody factors are too small to produce an observable 
effect in the airshower profile. Identical conclusions are reached by comparing the number of muons at various depths 
vs. the number of e + e~ pairs at the airshower maximum (right panel of Fig.[l|. 
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FIG. 1: Left panel: Number of e e~ pairs vs. slant depth for the longitudinal development of 50 ten-dimensional BH airshowers. 
The neutrino primary energy is E v = 10 19 TeV, the altitude of the first interaction depth is 10 km (slanth depth 780 g cm -2 ) 
and the zenith angle is 70°. BH airshowers without (with) QCD and emissivity effects are shown by black dashed (red dotted) 
curves. Right panel: Number of fi~ pairs at various atmospheric depths X m + AX vs. the number of e + e~ at the shower 
maximum. BH airshowers without (with) QCD and emissivity effects are shown by black filled (red empty) circles. The 
observation depth increases from left to right and top to bottom. 



III. ROTATING BLACK HOLES 



Up to now, simulations of BH airshowers have focused on Schwarzschild BHs @, [T^j. However, BHs created in 
collisions with nonzero impact parameter are expected to be spinning. Since the evaporation process depends on the 
BH angular momentum, airshowers initiated by spinning BHs could be significantly different from airshowers initiated 
by Schwarzschild BHs. If the graviton emissivity for rotating BHs is much higher than the emissivity of SM particles, 
only gravitons will be emitted, making the BH undetectable. Although particle emissivities for higher-dimensional 
rotating BHs are not fully known |l8|, the effects of rotation in the evaporation phase can be estimated from results 
in four dimensions and for higher-dimensional nonrotating BHs. 

In four dimensions, the evaporation phase of a spinning BH with large angular momentum is dominated by gravitons 
[l9| . The emissivity of spin-2 fields increases by a factor ~ 10 2 — 10 3 more than the emissivity of lower spin particles 
when the angular momentum increases from J = to the maximum value J m ax = M 2 , where M is the mass of the 
BH. For a random distribution of BH spins, the average increase in graviton emissivity is ~ 10 more than the other 
fields. Graviton emission also increases with the number of dimensions due to a higher number of spin-2 helicity states. 
This has been shown quantitatively in Ref. [TtJ for nonrotating BHs. The graviton-to-SM emission ratio increases 
from 1:10 3 in four dimensions to 1:4 in eleven dimensions. 

The results above suggest a larger graviton greybody factor for higher-dimensional, spinning BHs. The increase in 
graviton emissivity is specially relevant for ultra-spinning BHs. If most of the BHs produced in UHECR collisions are 
low-spinning, graviton emission is likely to increase on average by one order of magnitude more than the other fields. 
If most of the BHs are ultra-spinning, graviton emission could be enhanced by several orders of magnitude. It is thus 
crucial to determine the distribution of BH spins in airshowers. To this purpose, let us define the parameter a 0]: 

a = —MR' (3) 

where D is the number of dimensions and R is the BH radius. The radius of the BH is related to the mass and angular 
momentum by the relation (20| 



* = 7= 



1 + a 2 D-2 



8M T( 2 



Tl/CD-3) 

(4) 
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FIG. 2: Left Panel: Histogram of the number of events vs. a for the decay of a ten-dimensional rotating BH. Right Panel: 
Shower profiles for both rotating and nonrotating BHs. The black dashed (red dotted) curves denote rotating (nonrotating) 
BH events. 



We consider the model of Ref. [bj] for BH formation. The colliding particles are described by boosted Schwarzchild 
solutions at fixed energy (Aichelburg-Sexl shock waves) [2l[. The BH is formed when the two waves are superposed to 
form a trapped surface. The mass of the BH is related to the CM energy of the colliding particles, E cm , by M = E cm y, 
where y depends on the impact parameter b of the collision. The ratio J/M in Eq. is 

J _ xr 

where x — b/ro is the impact parameter normalized to ro = (47r£7 cm /r2.D_3) ^ and is the area of the unit 

sphere in D — 3 dimensions. The parameter a is the solution of the polynomial equation 

aP^ = C ±±*(*Y~* , (6) 



where 



C-M^r^-y CO 



The left panel of Fig. Ogives the distribution of a for 10,000 events in ten dimensions. Most of the BH are formed 
with small angular momentum. This is somehow expected because the BH production cross section is reduced by a 
factor of (1 + a 2 )~ 2 /( D ~ 3 ) compared to the nonrotating case JEHU- According to the discussion above, the increase 
in graviton emissivity relative to lower-spin fields can be estimated to be about one order of magnitude compared the 
nonrotating case. 

Simulations for spinning BHs with Tji 2 ~ IV 2 = 10 X (nonrotating Tn 2 1 IV2 ) show that most of the emission is in the 
form of gravitons. However, the number of observable secondaries for both rotating and nonrotating BH airshowers 
is stable and the difference in the profiles is statistically not significant (right panel of Fig. 2). In a typical event, the 
bulk of the collisional CM energy is carried by the nucleon remnant. Therefore, changes in particle emissivities have 
generally minor effects on the airshower characteristics. Rare events (< 10 %) are characterized by a very low energy 
output. This happens when the BH carries most of the collisional energy. In these cases, the nucleon remnant does 
not shower and the visible energy is highly reduced by the increased graviton emission from the BH. 
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IV. CONCLUSIONS 



This paper focused on two aspects of BH airshowers which had been neglected in previous studies: QCD and BH 
spin effects. The inclusion of these effects is important to check the stability of the airshower profiles and provide 
a more accurate template for observational searches. Event simulations based on the GROKE Monte Carlo show 
no change in the overall characteristics of the airshowers. Effects due to color conservation, initial- and final-state 
radiation, and different fragmentation models are washed out during the airshower development. BH spin effects in 
the airshower development are estimated to be typically small for two reasons: (i) most BHs are formed with low 
angular momentum and (ii) most of the CM energy is carried by the nucleon remnant. It should be stressed, however, 
that a conclusive statement on this issue requires the knowledge of the exact spinning BH greybody factors for all 
fields in higher-dimensions. 
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